function [discount, stipend] = bvp(lm)
discount=lm>=600;
stipend=(lm>=700)*800000;
end